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Abstract: We consider the recent proposal that the distribution of the difference between 
azimuthal angles of the two accompanying jets in gluon- fusion induced Higgs-plus- two- 
jet events at LHC reflects the CP of the Higgs boson produced. We point out that the 
hierarchy between the Higgs boson mass and the jet transverse energy makes this observable 
vulnerable to logarithmically enhanced higher-order perturbative corrections. We present 
an evolution equation that describes the scale variation of the azimuthal angular correlation 
for the two jets. The emission of extra partons leads to a significant suppression of the 
correlation. Using the HERWIG Monte Carlo event generator, we carry out a parton- 
shower analysis to confirm the findings. 
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1. Introduction 

Gluon fusion via the top quark loop provides the dominant mechanism for Higgs boson 
production at LHC. Recently, the analysis of the related (2 — > 3) cross section was carried 
out at the leading order in a>s(Mz) in ref. |IJ. 

The proposal is to tag two extra jets, the motivation being the elucidation of the 
differences/similarities between this process and the weak-boson-fusion process, which also 
comes accompanied with two extra jets. 

After imposing cuts on the jet momenta similar to the weak-boson-fusion selection 
criteria, the soft-gluon contribution is reduced and they obtain, as shown 1 in fig. [I], a 
striking correlation between the azimuthal angles <j) of the two jets, one of which is now 
in the forward direction while the other one is in the backward direction. Provided that 
the produced Higgs boson is CP-even, the distribution of the difference Acfi of the two 
azimuthal angles is peaked at A</> = 0, tt and falls to nearly zero at Acfi = ir/2. If the Higgs 
boson is CP-odd, although the result is not explicitly shown, the distribution is peaked at 
A(p = 7r/2 and falls to nearly zero at A(p = 0, it. 

In this work, we would like to point out one potential pitfall which seems to have been 
neglected in their study, namely that there are two scales in this problem. The higher scale 
is related to the Higgs boson production. For a leading order analysis we can set it to be 
the Higgs boson mass for convenience. The lower scale is related to the emission leading 
to the tagged jets and this is characterized by the jet transverse momenta. Because of the 
presence of two scales, the predictions of calculations at a finite perturbative order becomes 
sensitive to higher order corrections. 

We investigate this problem by first establishing an evolution equation that describes 
the scale variation of the azimuthal angular correlation coefficient. The large size of the rel- 
evant anomalous dimension implies that there is significant, up to one order of magnitude, 
reduction in the size of the correlation coefficient. 

x We thank the authors of ref. [EJ for their kind permission to reproduce the figure. 
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Figure 1: The distribution of the azimuthal angle between the two highest px jets, taken from 
ref. 0. Results shown are for the top-quark induced gluon-fusion process, with m t — 175 GeV and 
in the limit mt — > oo, and for the weak-boson- fusion process. 



Although the problem is formally due to the logarithm of the ratio of the two scales 
which enhance the higher order contributions, we find that the ratio needs not be so large 
for the effect of extra emission to become important. 

An alternative approach to investigating this problem is by a parton-shower level Monte 
Carlo simulation. We formulate this problem as a (2 — ► 1) cross section convoluted with 
the leading logarithmic parton shower using HERWIG [§], which includes the azimuthal 
spin correlations by default by using the algorithm of Collins and Knowles ||, ||, ||]. 

The results of the two approaches are in good agreement. 

We note that the weak-boson fusion mode is unaffected by extra emission, as the jet 
Pt scale is the only QCD scale present in this case. 

We present the evolution equation analysis in sect. [2] and its comparison with the 
parton-shower analysis in sect. ^. We present the conclusions in sect. |j. 

2. Evolution equation analysis 

Let us consider the evolution of a jet that is due to an initial state parton. In a Monte 
Carlo simulation, this is described by the backward evolution. 

At each stage of evolution, the probed parton has virtuality t, momentum fraction x, 
jet transverse momentum px according to some definition, and spin density p. We wish 
to measure the mean value of the component of p = pn that is aligned with the reference 
direction given by pr- The correlation arises in the first place because of the correlation, 
at the hard process level, between the planes of polarization of the gluons involved in the 
Higgs boson production. Hence the decorrelation between the plane of polarization of each 
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gluon and the direction of the related tagged jet is equivalent to the suppression of the 
azimuthal angular correlation between the two tagged jets. 

In the following, px and p are both two-component vectors. The spin density matrix 
Pij is given in terms of the vector p^ by pij = (1 + pk<rfj)/2. The third component of the 
vector p vanishes so long as the nucleon is unpolarized ||, |j. In terms of the distribution 
functions f(x,t,pr,p), we may write the scale variation as follows: 

d <p|| >{x,t) _ d f f{x,t,pr,p)p\\dp T dp 
dt dt J f(x,t,p T ,p)dp T dp 

1 f df(x,t,p T ,p) 



/ -Q t {P\\- <p\\>{x,t)}dp T dp. (2.1) 



The change in the distribution functions when the scale is raised from t to t + 5t is 
given by: 

5f(x,t,p T ,p) = £/m-<5/out, ( 2 - 2 ) 
5/out = j J dz^- ( ^-CF(z,(p,p)f(x,t,pT,p), (2.3) 

6 fin = ~T [ -^%p'CF(z,cj>,p')5(p- p(p>,z,cf>)) 
t J z 2tt 2tt 

x f(x/z,t,p' T (p T ,t,z,(j)),p'). (2.4) 

In eqn. <j> is the emission angle, z is the remaining energy fraction during the emission, 

such that the emitted energy fraction is 1 — z. p is determined by the following relation 
according to the notation of ref. j|] : 

p = (f 2 (z)n(0) + f 3 (z)p')/F(z, 0, (/), (2.5) 

where fi are various components of the splitting function shown in tab. [l], and n((j)) = 
(cos 2(f), sin 2(f)). In all of the above, C and F are the colour factor and the polarization 
dependent splitting function respectively, with F(z,<p,p') = P(z)/C + f ±{z)n{4>) • p 1 . 
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Table 1: The coefficients fi, taken from ref. [Q, modified to adopt the convention of ref. for the 
normalization. 



Combining the above with eqn. (|2.l] ), the contribution from the term Sf ou t vanishes 
and we obtain: 

d<P\\>(x,t) 1 fdf- m 

ai = Tom) J ~~dT ( p]r <n>^)) (2.6) 

1/t fdza s d(f) ^ h{z)n(<t>) + h(z)p' 



It I 

f(x,t) J z 2vr 2vr F(z,(f),p') 
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xf(x/z,t,p' T (p T ,t,z,(p),p') (p||- <p\\>(x,t)) 
_1 r (ha^d4 , , f(x/z,t,p' T ,p') 
t J z 2vr 2vr 9 Pt 



(2.7) 



/OM) 

[<p>(x, t)F(z, 0, p') - (/ 2 (*)n(0 + / 3 (z)^) ] • (cos 20 o , sin 20 o ). (2.8) 

In the last line, (j)Q is the azimuthal angle of pr- Between eqns. ( [2.7D and (|2.8j) , we have 
assumed that pt is defined in an additive manner, such that dpr = dp' T . Excepting this 
assumption, our discussion so far has been general. 

Now let us consider the case in which the direction of pr, or more generally the reference 
direction, is to a good approximation determined at one scale given, for instance, by the 
jet transverse momentum PTp an d we are interested in the soft emission that depolarizes 
the gluon between this scale and the hard process scale. 

In this case, 0o can be taken as constant. Integrating over cj>, by symmetry, the term 
proportional to n{4>) in eqn. ( fO| ) vanishes, as does the </> dependent term in F(z,cf),p'). 
We then have: 



d <p\\>(x,t) 



dint 



fdzas,, 
J ~^ dpdpT 



/ f(x/z,t,p' T ,p') 



f(x,t) 

P(z) <p\\>{x,t)-Ch{z)p\ 
dz a s f(x/z,t) 



z 2vr f(x,t) 
Hz) < P\\ > 0, t) - Cf 3 (z) < p|| > {x/z, t) 



(2.10) 



The pole at z — > 1 of P{z) cancels with the pole of fz{z). The expression of eqn. (|2.1C| ) 
is the difference between the change in the structure function f(x, t) given by the first term 
and the change in the spin density < pii > f(x, t) given by the second term. The two terms 
can be treated separately by introducing the plus prescription to account for the <5/ ou t term 
given by eqn. (^). We obtain: 



d <p\\>(x,t)f(x,t) 



dint 



dz as 
~2tt 



f(x/z,t) < Pl \>(x/z,t) Ch{z) 



df(x,t) 
dint 



(2.11) 
(2.12) 



Taking the x 3 moments of the above to obtain Mellin transforms, we convert the 
expressions into: 



d[< Pll >f](j,t) 
dint 

df(j,t) 



[< Pll > f](j,t) I dzJ-^CMx 



dint 

such that the anomalous dimensions are 

7<pn>/ 

7/ 



dzzi- 1 ^Cf 3 (z)+, 
dzzi- l ^-P(z), 



(2.13) 
(2.14) 

(2.15) 
(2.16) 
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respectively for the spin density and the parton distribution functions. 

As stated above, eqn. ( 2.1C| ) specifies the evolution of polarization in regions where 
extra emission does not alter the reference direction. For our purpose, the most interesting 
case is the soft/collinear gluon emission from the gluon line in between the hard process 
and the jet px scale. 

For small enough x, such that j ~ 1, the behaviour of the anomalous dimensions is 
controlled by the z — > region. If we choose the g — > gg splitting, the fact that P{z) has a 
pole at z — ► where as fa{z)+ does not, indicates that a large disparity arises in the two 
quantities. < pn > is the ratio of <p\\ > f and /, such that its running is governed by the 
difference of the above two anomalous dimensions. The relevant splitting function is: 



Ch{z) - P{z) 



-2C A 



(!-*)(! + *») 



(2.17) 



Exactly as in the case of the hadron multiplicity calculation || the expression of 
eqn. ( 2.16 ), when j = 1, is divergent as written, but a more careful analysis, taking into 
account the suppression of soft gluon emission due to coherence, leads to the result: 



A7O', a s ) 
[<p~^> f](l,t) 



Ca«5 j_ 
2tt 



1 



+ .... 



oc cxp 



exp 




(2.18) 
(2.19) 



Hence the evaluation of < pn > at the hard process scale ~ y/§, defining it to be the ratio 
of < pn > / and / individually integrated over the allowed, typically small, values of x, 
assuming perfect polarization at the lower scale ~ PTj, gives: 



<P||> (P 



Tj) 



cxp 



/ 6 Ptj I 6 s 



(2.20) 



Let us take the jet energy scale at around 20 GeV and the hard process scale at 100 
GeV. b = {11C a — 2nj)/127r at this order, and we can adopt for simplicity A ~ 200 MeV. 
We then obtain < pn >= 0.42. We are interested in the correlation between two jets, which 
is now given by the square of < p\\ >, and we obtain 0.18. Thus there is almost one order 
of magnitude dilution. 

There would be further dilution due to the modification of the px direction due to 
extra emission, and the imperfect polarization at the jet energy scale must also be taken 
into account. The parton-shower analysis, which we present in the following, incorporates 
all of these effects. 

3. Parton-shower simulation 



Our parton-shower simulation is based on the Monte Carlo event generator HERWIG 6.5 
Q. For the (2 — ► 1) hard process, 

99 - *, (3.1) 
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where $ is either a CP-even or a CP-odd Higgs boson with M$ = 120 GeV Q, we utilize 
the code for the MSSM H° and A production. We set tan/3 = 1, (3 - a = vr/2 and all 
supersymmetric masses equal to 10 TeV, such that the Higgs production is mostly due to 
the top quark loop only. The parton shower, including the azimuthal angular correlation 
|3|, [H, follows the HERWIG default, although we found two bugs which had to be corrected 
before the effect of the angular correlation could be seen. One bug affects the azimuthal 
angle of the branching q — ► gq, whereas the other, more serious, bug affects the Lorentz 
boost of the jets to the laboratory frame. The resultant code will be part of the next 
HERWIG sub-version. 

We have performed the analysis also incorporating the heavy quark initiated hard 
process QQ — ► but there was no significant alteration to the results presented herein. 

For jet reconstruction, we make use of the program GET JET ||. This uses a simplified 
version of the UAl jet algorithm with jet radius AR which we set equal to 0.4. 

For the Standard Model parameters, we adopted the default values in HERWIG. For 
the structure function, we have used the default set in HERWIG, namely the mean of the 
central gluon and higher gluon leading order structure functions of MRST98 || . We follow 
the numbers in ref. |l]] for the Higgs boson mass of 120 GeV and the jet selection criteria. 
We do not decay the Higgs boson and do not impose any cuts on its kinematics. We take 
the LHC centre-of-mass energy y/s = 14 TeV. 

The jet selection criteria, which we define to be applicable to the two jets with the 
highest transverse momenta PTj, are as follows: 

PTj > 20 GeV, \r]j\ < 5, AR 5j > 0.6 

\Arijj\ > 4.2, riji ■ Vj2 < 0, m j:j > 600 GeV. (3.2) 

Thus we are picking out events with one forward jet and one backward jet, with large 
combined invariant mass. 

The result of our parton-shower simulation on the azimuthal angular difference A(p, 
following the procedures outlined above, is shown in fig. |2[ The error-bars are for the statis- 
tics. We generated one million events in each case to compensate for the low acceptance. 
Since the total cross section is around 20 pb, this corresponds to about 50 fb -1 integrated 
luminosity. 

The difference in the normalization between the CP-even case and the CP-odd case is 
due to the different loop amplitudes depending on whether the produced Higgs boson is 
CP-even or CP-odd. The total cross section is similar to that shown in fig. [p. The points 
marked 'No correlation' corresponds to the same hard process as the CP-even case but with 
the azimuthal correlation turned off by means of the HERWIG option AZSPIN=. FALSE. 

The curves in fig. || are fitted by the Fourier series analysis using the expansion 
ao + aicos(2A^<). Our fit gives 101.2 + 13.0cos(2A</>) (fb) for the CP-even case and 
82.0 — 14.9 cos(2Ac/>) (fb) for the CP-odd case. With the azimuthal correlation switched off, 
we obtain 100.1 — 1.45 cos(2A</>) (fb) so that it is possible that the case without azimuthal 
correlation is not completely flat. We calculated the statistical error on the Fourier coeffi- 
cients to be about 1 fb for all coefficients. The coefficient in this case is found to be 1.24a 
away from being flat. 
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Figure 2: The distribution of the azimuthal angle between the two highest px jets, after imposing 
the cuts of eqn. (3.2). Simulation at the parton-shower level. The error-bars are for the statistics. 



The correlation is evidently diluted. As the distribution is almost proportional to 
1 + cos(2A</>) for the CP-even case in fig. [I], we may say that there is dilution by nearly 
one order of magnitude, in agreement with the analysis in sect. 

Next, let us look at the jet px scale dependence of the correlation coefficient. In order 
to account for the possibly intrinsically non-flat distribution, we calculate the mean of the 
CP-even and CP-odd correlation coefficients given by: 

1 

" 2 



a 



(3.3) 

CP-even \ "0 / CP-odd. 

For the evolution equation analysis, the correlation coefficient is calculated as the square 
of <p|| > given in eqn. ( 2.20| ). 

The result is shown in fig. ||, taking the jet px scale to be the minimum pxj for the 
HERWIG simulation and pxj in the evolution equation analysis as seen in eqn. ( 2.20 ). We 
rather arbitrarily chose A = 180 MeV in this case, equal to the value of the HERWIG 
variable QCDLAM. We set the higher scale ~ \/I = 120 GeV to correspond to the Higgs 
boson mass. 

We see that the HERWIG results are in good agreement with the evolution equation 
analysis in sect. ||. The running of the correlation coefficient at low jet energy scales is 
slower in HERWIG. This is natural when we consider the fact that the actual pTj is always 
higher than the imposed pxj cut, especially so when the cut is low. We verified this by 
taking the pTj > 10 GeV point and further imposing the constraint that neither of the 
tagged jets have pTj greater than 20 GeV. We obtain 0.000 ± 0.024 for the coefficient in 
this case. 

As mentioned in sect. there is further dilution in the HERWIG analysis compared 
to the evolution equation analysis due to the imperfect polarization at the jet px scale 
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Figure 3: Jet pt scale dependence of the correlation coefficient |ai/ao|, averaged over the CP-even 
and CP-odd cases. For the HERWIG numbers, the jet pr scale is taken to be the minimum Ptj> 
whereas for the evolution equation analysis, the jet px scale is as represented by pxj in eqn. (2.20). 
The error-bars on the HERWIG numbers are for the statistics. 



and the modification of the pr direction due to extra emission. The latter is expected to 
be important when the jet pt scale is low. The result in fig. || indicates that although 
these effects must be present, they do not significantly alter the behaviour derived by the 
evolution equation analysis. 
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Figure 4: The distribution of the azimuthal angle between the two highest px jets, after imposing 
the cuts of eqn. (3.2). Simulation at the hadron level, without soft underlying events. 



-8- 



Our analysis so far has been carried out at the parton-shower level. Let us now 
turn our attention to the effect of hadronization. In fig. f|, we show the numbers for the 
hadron level simulation corresponding to the cuts of eqn. ( |3.2| ), i.e., with pxj > 20 GeV. 
We have turned off soft underlying events. We have utilized the simplified calorimeter 
simulation in GET JET. For the electromagnetic calorimeter, we take the resolution to be 
a E /E = 10%/ yjE /GeV, and for the hadr onic calorimeter, we take the resolution to be 
oe/E = 50%/ \J E/GeV. We observe that except for the reduced acceptance rate, the 
result is similar to that shown in fig. |2[ The resulting coefficient < |ai/tto| > is consistent 
with the parton level simulation. 



4. Conclusions 



By means of an evolution equation analysis and the HERWIG parton shower, we have 
performed an all-order study of the jet azimuthal angular distribution of events containing 
a Higgs boson and at least two jets at LHC. 

The results of the two analyses are consistent with each other and predict that due to 
the emission of extra partons, the azimuthal angular correlation that arises in the fixed- 
order analysis is diluted by almost one order of magnitude. 

Although we have specialized to the case of Higgs boson production, our results are 
applicable to jets accompanying any hard process initiated by two gluons. 

We have verified that hadronization does not significantly modify our conclusions. 
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